Dynamically generating arbitrary spin-orbit couplings for neutral atoms 
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Spin-orbit coupling (SOC) can give rise to interesting physics, from spin Hall to topological 
insulators, normally in condensed matter systems. Recently, this topical area has extended into 
atomic quantum gases in searching for artificial/synthetic gauge potentials. The prospects of tunable 
interaction and quantum state control promote neutral atoms as nature's quantum emulators for 
SOC. Y.-J. Lin et al. recently demonstrated a special form of the SOC k x a y : which they interpret 
as an equal superposition of Rashba and Dresselhaus couplings, in bose condensed atoms [Nature 
(London) 471, 83 (2011)]. This work reports an idea capable of implementing arbitrary forms of 
SOC by switching between two pairs of Raman laser pulses like that used by Lin et al. While one 
pair affects k x a y for some time, a second pair creates k y a y over other times with Raman pulses 
from different directions and a subsequent spin rotation into ±k y a x . With sufficient many pulses, 
the effective actions from different durations are small and accumulate in the same exponent despite 
that k x a y and ±k y a x do not commute. Our scheme involves no added complication, and can be 
demonstrated within current experiments. It applies equally to bosonic or fermionic atoms. 

PACS numbers: 03.75.Mn, 67.85.Fg, 67.85. Jk 



Introduction. Atomic quantum gases are increasingly 
viewed as favored model systems for emulating condensed 
matter physics. Optical lattices resulting from ac Stack 
shifts to atomic levels, are easily implemented with co- 
herent laser beams, which confine atoms like electrons in 
solid states. An interesting topic concerns strong correla- 
tions as in integer/fractional quantum Hall effect and the 
analogous spin Hall effect. The standard description for 
the former involves U(l) Abelian gauge fields, which can 
be simulated in neutral atoms through rotation 0, 3 or 
adiabatic translations in far-off- resonant laser fields 3|— 16|] . 
Non- Abelian gauge fields, e.g., as in spin-orbit coupling 
(SOC) [7- 11 1, enable richer possibilities like fractional 
quantum Hall states. As a result, active researches are 
targeting the implementations of (SOC) in simple atomic 
systems. 

For atoms with multiple internal states, or (pseudo- 
spin) spinor degrees of freedom, SOC changes single par- 
ticle spectra and competes with density-density or spin- 
dependent interactions, (i.e., spin-exchange and singlet- 
pairing interactions). Str ong correlations often lead to 



exotic ground states [12| - |21| . such as the plane- wave 



phase and the striped phase discovered recently in pscudo 
spin- 1/2 (H-EDl or spin-1 condensates [l3j]. Other exam- 
ples offer the triangular-latticed phase or square-latticed 
phase in spin-2 condensates with axisymmetric SOC 
[lil Il7j|. In a recent experiment, the JQI group of Spiel- 
man observed both Abelian @ and non- Abelian [§| gauge 
fields in a pseudo spin-1/2 atomic Bose gas, albeit in 
a special form oc k x a y of SOC, which is an equally 
weighted sum of Rashba (oc k x o~ y — k y o~ x ) and Dressel- 
haus (oc k x a y + k y a x ) couplings More generally, a 
SOC form of continuous rotation symmetry, or an arbi- 
trary weighted sum of Rashba and Dresselhaus couplings, 
exists in solid-state materials. 



Several existing theoretical proposals are capable of 
implementing SOC with rotation symmetry in laser atom 
coupled models. For instance, in a tripod scheme 0, 
when one-photon resonant couplings between the three 
lower-energy states and a higher-energy one are allowed, 
two dark states emerge, although spontaneous emission is 
always a cause of concern in this case. D. L. Campbell et 
al. [1(3 ] proposed an alternative scheme by cyclically cou- 
pling three or four ground or metastable internal states. 
With sufficient laser intensities, the above induced SOCs 
can possess a continuous rotation symmetry. Another 
scheme by Jay D. Sau et al. ll| employs an effective 
two-dimensional periodic potential created from two laser 
beams and their reflected lights propagating along x and 
y directions in 40 K atoms. In the limit of small Raman 
coupling, their corresponding effective SOC is of a pure 
Rashba type in the first Brillouin zone. 

In this Letter, we describe a dynamic approach for im- 
plementing rotational symmetric SOC of arbitrary forms 
within a pseudo-spin 1/2 atomic system. We adopt the 
JQI model and start from the simple SOC they proposed 
and recently demonstrated The key to our idea is 
optimal control theory applied with repeated laser pulses 
to rotate atomic pseudo-spins. Our idea works for both 
atomic fermions and bosons, and can be easily adopted 
to other atomic models. Thus it constitutes a power- 
ful new direction for engineering synthetic atomic gauge 
potentials. 

The equally weighted sum of Rashba and Dresselhaus 



types SOC of k x a y 



can be rotated into a form 



±fcj,(X 



by performing single atom spin rotation through 
a Rabi pulse. Such a coherent control idea when re- 
peated over time, can realize k x o~ y and ±k y o~ x types SOC 
in subsequent time intervals of duration St. The result- 
ing dynamics is then described respectively by an effec- 
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tivc Hamiltonian with pure Rashba or Drcssclhaus SOC 
with the first order approximation for small St. The ac- 
companied change of atomic momentum, can be nullified 
through a variety of means as we describe below step by 
step. We start with a review of the experiment by Y.-J. 
Lin et al. Q, which helps to introduce our idea. 

The JQI protocol. Consider a F = 1 atomic Bose- 
Einstein condensate (BEC) under a bias magnetic field 
along z located at the intersection of two Raman laser 
beams propagating along y + z and — y + z, with an- 
gular frequencies lol and u>l + Awl, respectively. The 



two laser beams affect two photon resonant Raman cou- 
pling (flu) between nearby ground Zeeman states, far 
detuned from the excited states. Effectively, such a 
coupling scheme produces an artificial magnetic field 
along the cc-axis direction with the resulting Hamiltonian 
Hr = VtRF x cos{2kLy + Aw^i), where F X)V)Z are 3x3 
spin-1 matrices, &l = ^/2n/X with A is the laser wave- 
length, and El = Ti 2 k\/2m, the unit of photon recoil 
energy. In explicit forms, after adiabatically eliminating 
excited states, the total Hamiltonian becomes 



2£2 



H* = 



2m 



E + 
E 
E_ 



V2 



cos(2fc L y + Auj L t) 

cos(2fciy + Auj^t) cos(2kLy + Aw^t) 



(1) 







cos(2fci,y + AwLt) 







where E+, Eq and £L are Zeeman (eigen-) energies of 
Mf = 1, 0, —I spin states, respectively. Under the rotat- 
ing wave approximation, it turns into 
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Further introduce a frame transformation: 
g-iF^AuLi^ where "0 and fp are the wave functions in 
the laboratory and transformed frames, respectively, we 
arrive at the Hamiltonian 



Ho 



h 2 k 2 

2m 



2Huj, 



, 3(5/2 
5/2 





F x cos(2k L y) 



F y sm(2k L y), 



8/2 
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where Hluz = E— — Eq, HAlol = hu>z + 5, Eq — E + = 
hujz — 2huj q , 5 is detuning and hui q is the quadratic Zee- 
man shift. When hoj q is sufficiently large and the Ra- 
man coupling SI = Qr/\/2 is small, we neglect the state 
\Mp = 1) and a constant term Eq — 5/2. The effec- 
tive Hamiltonian for the remaining two nearly degenerate 
states becomes 
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where the second line shows an explicit SOC term oc k y a z 
when viewed after a unitary transformation. 



(a) 



(b) 




FIG. 1: (Color online), (a) A schematic illustration of the 
JQI implementation for SOC [jj, where a F = 1 atomic Bose- 
Einstein condensate interacts with a bias magnetic field along 
z and two Raman laser beams propagating along y + z and 
— y + z, with angular frequencies lol and ujl + Aljl, respec- 
tively, (b) LEFT: Linear Zeeman shifts of the three hyper- 
fine spin states. MIDDLE: Zeeman shifts of the three hy- 
perfine spin states including both linear and quadratic terms. 
RIGHT: Zeeman shifts in the rotating frame (with frequency 
AlulFz) of the pseudo-spin pointing along z . 

Dynamically generating arbitrary SOC. Our protocol 
for implementing the Rashba type SOC is illustrated be- 
low in Fig. [5] It relies on our ability of being able to 
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switch atomic pseudo-spin from along z- to along y-axis 
(and vice versa) using Raman pulses. In the first half pe- 
riod, Raman lasers L\ and L 2 are turned on. In the sec- 
ond half, L 3 and L 2 are turned on instead. L 3 is the same 
as L\ except it propagates along opposite direction. At 
the middle point, we pulse on an extra ir/2 pulse to rotate 
the pseudo-spin from y- to z-axis, described by the oper- 
ator exp[— i(cr x /2)iv/2] in the transformed frame, or the 
operator exp{iF z Au)Lt] exp[— i{a x /2)-K /2] exp[— iF z Auj^t] 
in the lab frame; in the end of each period, we pulse on an 
— 7r/2 pulse for the reverse rotation. Both spin rotation 
pulses can be accomplished with cither Raman coupling 
from appropriately detuned lasers or rf plus microwave 
coupling between the two remaining internal states. 

In the first half, the system is then governed by H3 of 
the Eq. © as in the JQI experiment |9j. In the second 
half, the Hamiltonian in the transformed frame following 
that in the Eq. ([3]) becomes 

h 2-r 2 / 2hu q + 35/2 
H' B = — — + 5/2 I + E - 5/2 

Zm \ -5/2 
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which completes one period of our prescribed protocol. 
For large TiLo q and small ft, the same condition as in Rcf. 

the effective Hamiltonian for the reduced two-state 
model becomes 
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The pair of 7r/2 pulse (before) and —ir/2 pulse (after) 
affects a unitary transformation 

e i( CTx /2)7r/2£/ e -i( CTai /2),r/2 

Kb 2 5 ft . , ft . , M 
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to compensate for photon recoils 



±— pulse to rotate atomic spin 



FIG. 2: (Color online), (a) A schematic diagram for dynam- 
ically generating arbitrary spin-orbit coupling, (b) The pulse 
sequence used to implement Rashba type SOC. The blue and 
cyan ones are suitable momentum impulses used to compen- 
sate for the unwanted photon recoils, which can be accom- 
plished with artificial or real inhomogeneous magnetic fields 
or suitably arranged state dependent Bragg pulses. The red 
±7r/2 pulse rotates atomic spin. 
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direction and with a spatial gradient. Thus they can 
be nullified by real magnetic field gradients or synthetic 
magnetic field gradients generated from spatial depen- 
dent ac Stark shifts. For instance, U\ is compensated for 
by a magnetic field pointing along z-axis and a spatial 
gradient (B') along y-axis, with an adjustable impulse 
over St' where E' = ~[iB' is the appropriate Zeeman en- 
ergy gradient. After the control pulse St, the sign of B' 
is changed to affect a second impulse, which then leads 
to the following 
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In suitably transformed frames, respectively with U± = 
e -ik L ycr z and U2 = e ik L 2<r v ^ Eqs gj) and q r eveal explicit 

SOC terms k y a z and k z a y . They cannot, however, be 
simply added together in the forms above. To combine 
the above two halves into a single Rashba or Dresselhaus 
type SOC, we have to eliminate these unitary transfor- 
mations. Both XJ\ and U 2 corresponds to spin dependent 
phase shifts, they can be viewed as from the impulse of 
an artificial or real small magnetic field along a suitable 
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provided E' St' = 2hki Jl where we assume E' is strong 
enough so that we can neglect the contribution from H3 
during the short pulse St' (<§; St) . The effective two-state 
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dynamics is then approximately govern by 
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apart from a overall phase term involving a constant en- 
ergy in the exponent. Similarly, U<i is nullified as well, 
resulting in 
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For the special case of Rashba SOC, the suggested 
pulse sequence are illustrated in Fig. [U(b), where the blue 
and cyan ones are suitable momentum impulses for com- 
pensating the unwanted momentum recoils in the first 
and second half cycles respectively. The red pairs are 
±7r/2 pulses for rotating the pseudo-spin. If the it/ 2 one 
precedes the — 7r/2 pulse, we find in one period T = 2St, 
the total evolution operator under two-state approxima- 
tion is given by 
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According to the Floquct theorem, the quasienergy e 
of time-periodic system is derived from dot [L/(T, 0) — 
e~ teT ] = 0. Then from Eq. (fT2|) . we can easily infer that 
under first order of T approximation, the quasienergy 
of our system is the same as the spectra of that with 
Rashba SOC. Reversing the two red ±7r/2 pulses intro- 
duces a minus sign " — " , the Rashba SOC then changes 
into Drcssclhaus SOC. By adjusting the timing constant 
St, we can extend the above discussion to SOC of arbi- 
trary form f3(k y cr z - k z a y ) + y/\ - |/3| 2 (k y a z + k z a y ). 
The steady state of the effective system Hamiltonian is 
reached due to elastic atomic collisions. Although in the 
simplest case, one period of the control protocol is often 
sufficient, the actual implementation can aim at a higher 
precision of the effective SOC Hamiltonian by increasing 
the number of cycles, or simple reducing St. 

A magnetic field gradient was first used in Ref. [6( for 
implementing Abclian gauge fields with neutral atoms. 
However, since real static B-field is subjected to the 
Maxwell's equations V ■ B = and VxB = 0, one cannot 
simply obtain a linear gradient along one direction, e.g. 
a B-ficld like B = (Bq — by)e y is illegitimate because of 
its non-vanishing divergence. The simplest linear gradi- 
ent B-field, therefore needs to have two components, like 
that of the commonly used two dimensional quadruple 
field, B = bxe x — bye y . When the system is of a reduced 
dimension, not including x-direction as in Ref. 0, one 
is then equipped with a one-dimensional B-ficld gradient 



B = —bye y , which is equivalent to B = —bye z upon an 
axis rotation, A from precisely needed for implementing 
U\ above. 

Likewise, the above B-field gradient can be simu- 
lated using ac stark shifts from position dependent laser 
fields far off resonant coupled to the two states form- 
ing an atomic pseudo-spin. Assuming a one-photon 
resonant coupling Rabi frequency Qj,(f) and a detun- 
ing Al = f^L — Qq, the ac Stark shift takes the form 
oc cr z |f2| i (r^|/(2Ai) oc o z lh{r)l 'Al- A linear spatial gra- 
dient can thus be affected with a laser intensity gradi- 
ent, which can be implemented using many methods, 
including the use of a gradient neutral density filter. 
Stronger gradients arise from interfering several waves 
forming a standing wave, e.g., with Ihif) ~ cos 2 (qy/2) = 
(1 + cosqy)/2, linear gradient is oc ±qy around the nodal 
points of cos qy. 

More generally, the state dependent gradients can be 
engineered to couple states in the same Zeeman manifold. 
For example, the above ac Stark shifts from one-photon 
coupling can be substituted with two-photon Raman cou- 
pling with suitable differential detuning, like in Bragg 
scattering, which then implements impulses oc ±yer 2 , 
oc ±za y or oc ±za z . 

Summarizing We present a coherent control protocol 
capable of realizing the Rashba type SOC in a pseudo- 
spin 1/2 atomic quantum gas [9j. For most systems, 
our protocol can be implemented in one cycle, involving 



5 



two separate resonant Raman coupling. More elaborate 
forms are possible with multiple control pulses. When 
more than one control cycle is implemented, we can fur- 
ther enhance the precision and strength of the SOC, or 
the corresponding artificially created gauge potentials. In 
addition, the scheme we suggest is independent of quan- 
tum statistics of atoms, thus can be adopted to fermionic 
atoms as well. Our idea thus opens the door for dy- 
namically implementing artificial gauge potentials in cold 
atomic systems based on coherent control theory. 

Finally, we compare our idea with two previous 

schemes [1(1 EH . In Ref. , three and four laser fields 
are needed, cyclically coupled to three or four internal 
states. Nearly pure Rashba or Dresselhaus SOC then 
results respectively in the limit of large intensity laser 
fields. It remains open to find a suitable experimental 
system. In Ref. , along each axis of x- and y- two 
lasers with different frequency and their respective re- 
flections are needed. Only in the far-detuned and small 
Raman coupling limit, Rashba or Dresselhaus SOC 
can be implemented, which results in a relatively small 
SOC, proportional to Qr. Our scheme, however, takes 
the full advantage of the JQI protocol [§]. By simply 
turning on several pulses, and making use of the coher- 
ent control, we can dynamically generate arbitrary SOC 
for neutral atoms. 
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